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We study the spontaneous dissociation of diatomic molecules produced in cold atomic gases via 
magnetically tunable Feshbach resonances. We provide a universal formula for the lifetime of these 
molecules that relates their decay to the scattering length and the loss rate constant for inelastic 
spin relaxation. Our universal treatment as well as our exact coupled channels calculations for *^Rb 
dimers predict a suppression of the decay over several orders of magnitude when the scattering 
length is increased. Our predictions are in good agreement with recent measurements of the lifetime 

of «5Rb2. 
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The discovery of Feshbach resonances in cold gases pro- 
vides the unique opportunity to widely tune the inter- 
atomic interactions via magnetic fields near a singularity 
of the scattering length. This singularity is due to the 
near degeneracy of the very small collision energy of the 
cold atoms with the binding energy of an extremely loose 
diatomic molecular state for magnetic fields at which the 
scattering length is positive. A number of recent exper- 
iments, as summarised in the accompanying experimen- 
tal paper of Thompson et al. , have taken advantage 
of this virtual energy match to populate this molecular 
state. The interacting pairs of atoms in the first stud- 
ies of these cold Feshbach molecules 0, 1^ were subject 
to deeply inelastic spin relaxation collisions involving the 
dcexcitation of at least one of the colliding atoms. Spin 
relaxation can affect both free Q and bound atom pairs. 

In this letter wc predict the spontaneous dissociation of 
Feshbach molecules due to spin relaxation. We provide 
a universal treatment of the molecular decay by which 
we relate the molecular lifetime r to the spatial extent of 
the bound state wave function and the loss rate constant 
for inelastic spin relaxation collisions. Both the size of 
the molecule and the loss rate constant exhibit a pro- 
nounced resonance enhancement at the singularity of the 
scattering length. General considerations on the nature 
of this enhancement allow us to predict the functional 
form of the magnetic field dependence of r. We show 
that the near resonant molecular lifetime increases with 
increasing scattering lengths and can therefore be varied 
over several orders of magnitude in related experiments. 
We demonstrate the predictive power of our universal 
treatment in a comparison with exact coupled channels 
calculations of the lifetime of the ^^Rb dimers produced 
in the experiments of Refs. The exact calcula- 

tions show that T varies in the remarkable range between 
only 100 fis and roughly 30 ms in the experimental range 
of magnetic field strengths. Our studies indicate that 
it is the large spatial extent of the near resonant bound 



state wave function that can provide the crucial stability 
of the Feshbach molecules with respect to their sponta- 
neous decay. 

In the following, we shall consider pairs of atoms in 
a cold gas that interact via s waves. We denote the bi- 
nary scattering channel of a pair of asymptotically free 
atoms as the entrance channel and choose the zero of en- 
ergy at its dissociation threshold. This threshold is de- 
termined by the Zeeman energy of the separated atoms. 
The Hamiltonian H of the relative motion can then be 
divided into the contribution Hd of the closed channels 
with dissociation thresholds at or above zero energy (in- 
cluding the entrance channel), the contribution Hop of 
the open channels with negative threshold energies, as 
well as the weak coupling Hi^t between the closed and 
open channels: 
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We shall consider magnetic field strengths on the side 
of positive scattering lengths of a zero energy resonance 
(i.e. a singularity of the scattering length). The large pos- 
itive values of the scattering length imply the existence of 
a weakly bound (metastable) molecular state 4>'^i, which 
corresponds to the highest excited vibrational {v = — 1) 
bound state of Hc\. The wave function thus fulfils 
the stationary Schrodinger equation Hc\(j)'^i = E^i(f)'^i, 
where E'^^ is the negative binding energy. Fermi's golden 
rule determines the exponential decay of the Feshbach 
molecule by a rate 
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Here Gop{E^^i+iO) = [E"}^ + iQ - Hop) ^ is the Green's 
function of the open binary scattering channels, whose 
energy argument E'^i + iO indicates that E'^i is ap- 
proached from the upper half of the complex plane. As 
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\E'^-^ I is much smaller than typical transition energies, it 
can be neglected in the argument of the Green's function. 

Fermi's golden rule also determines the loss rate con- 
stant for inelastic spin relaxation collisions to be 
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Here 
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is the zero energy scattering state associ- 



ated with the Hamiltonian Hd, which fulfils the sta- 
tionary Schrodinger equation Hci4>l)^ci ~ 0. The prcf- 
actors in Eq. JSJl correspond to a cold thermal gas of 
identical Bosons. The associated rate equation reads 0| 
N{t) = -K2{n{t))N{t), where N{t) is the number of 
atoms and {n{t)) is their average density. 

At magnetic field strengths in the vicinity of a zero 
energy resonance, 4'o'c\ ^^"^ ^'"^ determined by their 
wave functions in the entrance channel and in a closed 
channel strongly coupled to it 0, • This strong inter- 
channel coupling is due to the near degeneracy of the 
energy E^cs{B) of a closed channel vibrational state (the 
bare Feshbach resonance level) (jj^cs with the dissociation 
threshold of the entrance channel. The restricted two- 
body Hamiltonian Hd of the entrance channel and the 
closed channel can be effectively described by: 
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Here r denotes the relative distance between the atoms 
and m is twice their reduced mass, B is the magnetic 
field strength, Vbg{r) is the background scattering po- 
tential of the entrance channel, and W{r) provides the 
inter-channel coupling. The closed channel potential 
Vc\{B,r) supports the resonance state, i.e. [— ?i^V^/m + 
Fci(B,r)]0,.cs(r) = E,.cs{B)(f>rcs{r)- The dissociation 
threshold of Vc\{B,r) is determined by the energy of a 
pair of non-interacting atoms in the closed channel. The 
relative Zceman energy shift between the two channels 
as well as the bare energy E,:cs{B) can be tuned by vary- 
ing the magnetic field strength B. We denote by B^cs 
the magnetic field strength, at which E,-cs{B) crosses the 
dissociation threshold energy of the entrance channel, 
i.e. i?ros(-Bros) = 0. The bare energy E^csiB) varies vir- 
tually linearly in B, and an expansion about B^es yields 

/ijcs is the difference in 
magnetic moment between the Feshbach resonance state 
and a pair of atoms in the entrance channel. 

Under the assumption that the spatial configuration 
of a pair of atoms in the closed channel is restricted to 
the resonance state (/)res ('''), the dressed highest excited 
vibrational bound state ip'^^ of Hd is given by 0: 
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Here Gbg{E^^^) = iE^_}i + h^V^ /m-Vbg)-^ is the Green's 
function associated with the entrance channel, and the 



factor 1/Af assures that i/)^]^ is unit normalised. The as- 
sociated binding energy E'^i is determined by E'^i = 
E,osiB) + {<j)res\WGbgiE'_}^)W\^res)- The stroug coupling 
between the entrance channel and the bare Feshbach res- 
onance state (pres {r) also determines the zero energy scat- 
tering state 



j)^d of -^ci by the general formula 0: 
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Here 



uq denotes the bare zero energy state associ- 
ated with the background scattering, i.e. [—h'^\7'^/m + 
Vbg(r)]4+Jg(r) = 0, and the amplitude A{B) in Eq. ® 
is given by 
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The scattering length a{B) can then be obtained from 
the asymptotic behaviour (27r?i)^^/^ [1 — a{B)/r] of the 
entrance channel component of Eq. © at large inter- 
atomic distances r. This yields 
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where AB = {2nh)^\{(f>,^s\W\(f>l^^^)\^m/iA7Th^abgtires) is 
the resonance width and Cbg is the background scattering 
length. The singularity ofA{B) thus determines the mea- 
surable position Bq of the zero energy resonance, which 
is distinct from the zero B-^cs of the bare energy Ej-cs{B) 
by the shift Bq - B^cs = -((?!>rcs|M^Gbg(0)M^|(^rcs)/Mrcs of 
the denominator on the right hand side of Eq. . 

At near resonant magnetic field strengths the binding 
energy is determined by E'^^ = —h^/ima^) (cf., e.g., 
Refs. H, 0)- The energy argument of the Green's 
function in Eq. ((SJ, therefore, vanishes at large positive 
scattering lengths. In this limit the background scatter- 
ing contribution 0o^g to the right hand side of Eq. © 
also becomes negligible for inter-atomic distances on the 
length scale set by the range of the inter-channel coupling 
ffint in Eqs. ^ and (0). Equations (jSj and Q thus co- 
incide up to a factor MA{B). and Eqs. (|2J) and (j^l yield 
K2iB)/j{B) = 2{2Trhf\A{B)\^Af^iB). The near reso- 
nant amplitude A(B) is readily obtained from the reso- 
nance width and shift of Eq. © , while the normalisation 
constant Af'^{B) ~ {^^f,sAB)ma\,ga{B) / {2h'^) can be de- 
termined from {(tfli\(tfli) = 10. The molecular lifetime 
t{B) — 1/7(5) is then given by the universal formula 

t{B) = iTTa^{B)/K2{B) = [4/A'2(B)]4^(r3)/3. (9) 



(5) Here (r^) is the expectation value of for the near reso- 



nant wave function (l)-i{r) = exp(— r/a)/(r-\/27ra) 0|1|- 
General considerations [9j show that in the presence of 
open decay channels the resonance enhanced zero energy 
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binary elastic and inelastic collision cross sections can be 
described in terms of a complex scattering length a{B) — 
ib{B), whose imaginary part is related to the loss rate 
constant by 
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K2{B) = 16Trnb{B)/r 
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Its real part is well approximated by Eq. ijHJ), except that 
a?{B) approaches a large but finite value at B = Bq, 
which is determined by the decay width. We shall pre- 
suppose in the following that only the bare resonance 
state (j)T:os decays to a set of final channels with a total 
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is the associated lifetime. 



The imaginary part of the complex scattering length is 
then given by a Lorentzian function 
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Inserting Eqs. (jHJ, COl and lfTT|) into Eq. © then yields 



t(B) = r,. 



' AT? ^ 



r,.es/(B), (12) 



where we have introduced x = {B — Bo)/AB. This for- 
mula provides the general functional form of t{B) as well 
as its order of magnitude in terms of the constant factor 
Tics, which can be obtained either from a measurement 
of K2{B), at a single magnetic field strength, or from 
coupled channels calculations. 

We shall demonstrate the predictive power of Eqs. © 
and (|12|l for the example of the Feshbach molecules as- 
sociated recently from pairs of ^^Rb atoms in a Bose- 
Einstein condensate P, 0| at magnetic field strengths 
of about 15.5 mT. The nuclear spin / = 5/2 of ^^Rb 
gives rise to two energetically relevant hyperfine levels 
with total spin quantum numbers / = 2 and 3. A ho- 
mogeneous magnetic field B splits the Zeeman sublevels 
into a non-degenerate set of levels {f,mf), labelled by 
the projection quantum number m/ along the axis of 
the magnetic field and the / value with which it cor- 
relates adiabatically at zero field. In the experiments 
0, 01 the atoms were prepared in the (/ = 2, m/ = —2) 
internal state. The binary asymptotic scattering chan- 
nels are then characterised by pairs of internal quantum 
numbers (/, m/; /', m^), as well as the orbital angular 
momentum quantum number i associated with the rel- 
ative motion of the atom pair and its projection quan- 
tum number mi. The s wave {£ = 0, = 0) entrance 
channel is thus characterised just by the internal quan- 
tum numbers (2,— 2; 2,— 2). The cylindrical symmetry 
of the setup implies that the projection quantum num- 
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' iTif = —4 of the total angular momentum 



is conserved in a binary collision. 

We have applied the coupled channels approach of 
Ref. to determine the scattering length a{B) (see 
Fig. GJ, including spin exchange interactions, which cou- 
ple the entrance channel to the four closed s wave scatter- 
ing channels (3, -3; 2,-1), (3, -2; 2, -2), (3,-1; 3, -3) 
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FIG. 1; Magnetic field dependence of the real (solid curve) 
and imaginary (dashed curve) parts of the complex scattering 
length a{B)—ib{B). The dotted line indicates the length scale 
Hi rvdw = {mC^lh^f^^ 12 set by the long range asymptotic 
—C^jr^ van der Waals tail of the background scattering po- 
tential. We note that h{S) is multiplied by a factor of 1000. 



and (3,— 2; 3,— 2). A fit to Eq. © predicts the param- 
eters So = 15.520 mT, AS = 1.065 mT, and Obg = 
-484.1 ao (ao = 0.0529177 nm) for the experimentally 
well characterised 15.5 mT zero energy resonance of 
^^Rb 0. The decay of the Feshbach molecules is due 
to spin-dipole interactions, which couple s waves to d 
waves. Among the 23 d wave channels coupled to the en- 
trance channel only three are open. These d wave decay 
exit channels are characterised by the internal quantum 
numbers (2, -2; 2,-1), (2, -2; 2,0) and (2,-1; 2,-1) in 
order of increasing energy release. The transition en- 
ergies are all on the order of several mK (in units of 
/CB = 1.3806505 X lO"^^ j/k), which implies that the 
decay products are lost from an atom trap. Including 
spin-dipolc interactions in the coupled channels calcula- 
tions shifts the predicted position Bq by less than 1 /iT 
but it provides the small imaginary part h{B) of the scat- 
tering length (see Fig.^l, which determines the loss rate 
constant K2{B) by Eq. (P|l . 

Given the magnetic moment difference /iros/^ = 34.66 
MHz/mT, a fit of Eq. ^ to the curve of h(B) in Fig. [2 
determines the lifetime of the bare resonance level to be 
I/710S — 32 /is. The solid curve in Fig. |5] shows 



the associated universal estimates of the lifetimes t(B) of 
the dressed ®^Rb2 Feshbach molecules 0, 01 determined 
from Eq. JHl), using the exact loss rate constant K2{B) 
[cf. Eq. H1U|I and Fig.^. The dashed curve indicates the 
predictions of Eq. (|12() . We expect these estimates to 
be accurate between Bq and about 16.0 mT, where the 
scattering length a{B) by far exceeds the length scale 
Tvdw = 84 ao set by the long range asymptotic behaviour 
of the van der Waals interaction (see Fig.^. This trend is 
confirmed by our direct calculations of the molecular life- 
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FIG. 2: Two-body decay lifetime r(B) of the **^Rb2 Feshbach 
molecules 0, I3 versus magnetic field strength B near 15.5 
mT. The squares indicate the experimental data of Thompson 
et al. The circles indicate exact predictions of coupled 

channels calculations, whereas the solid and dotted curves 
show the results of the universal formulae (|5J and 1121 . re- 
spectively. The magnitude of the bare state lifetime Tros is 
indicated by the dotted line. For the purpose of compari- 
son, we have slightly corrected the axis of the magnetic field 
strength in such a way that the theoretical resonance position 
recovers the measured value of Bo = 15.5041 mT 1^. 

time shown in Fig.[5| The exact lifetime can be obtained 
from the energy dependence of the inelastic scattering 
cross section between any pair of open decay channels. 
We have chosen transitions between the open d wave 
channels associated with the internal atomic quantum 
numbers (2, -2; 2, -1) and (2, -2; 2, 0). Fitting the stan- 
dard Breit-Wigner formula to the resonance peak at the 
binding energy E\,{B) 0, of the metastable Feshbach 
molecules then determines the decay rate ^{B) = \/t{B) 
through its width. 

The squares in Fig. indicate lifetime measurements 
by Thompson et al. [ij in which the dimer molecules 
were produced in a dilute vapour of ^^Rb atoms at a tem- 
perature of 30 nK using a Feshbach resonance crossing 
technique (see, e.g., Refs. 0,0)- Our theoretical predic- 
tions agree with the experimental lifetimes in both their 
systematic trends and magnitudes. Figure [21 shows the 
large differences in the order of magnitude of t{B) from 
roughly 30 ms to only 100 /xs in the experimentally ac- 
cessible 0] range of magnetic field strengths from about 
15.6 to 16.2 mT. The lifetime correlates with the admix- 
ture of the bare resonance level to the dressed molecular 
state, which, in turn, is related to the spatial extent of 
the molecules . Equation @ can indeed be interpreted 
in the intuitive way that 47r(r'^)/3 sets the scale of the 
volume confining the bound atom pair, while K2{B)/'i 



is the event rate constant for molecular decay (which is 
the same as that for the inelastic collision of a pair of 
identical Bose condensed atoms 0). The large extent 
of the molecular wave function at near resonant mag- 
netic field strengths thus protects the Feshbach molecules 
against their spontaneous decay. We caution, however, 
that purely two-body theory may eventually become in- 
adequate sufficiently close to resonance; for example, the 
mean distance between atoms at the experimental Q 
peak density of 6.6 x 10^^ cm~^ is only 2 x 10* oq, which is 
comparable to the estimated bond length of the molecules 
111 of about 3 X 10^ oq at 15.6 mT. 

Our general theory applies to a variety of Feshbach 
molecules with open decay exit channels Il2l . Among 
these species are the *°K2 and ^Li2 |l3l | molecules 
produced in Fermi gases in the s wave entrance channels 
associated with the pairs of atomic quantum numbers 
(9/2, -9/2; 9/2, -5/2) and (1/2, 1/2; 3/2, -3/2), respec- 
tively. Measurements of two-body decay lifetimes could 
allow the characterisation of the spin-dipole interaction 
between various atomic species with an unprecedented 
precision. Our universal estimate determines the con- 
ditions for the stability of Feshbach molecules, which is 
crucial to all future studies of cold molecular gases in the 
presence of open decay exit channels. 
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